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Abstract 



This work considers special types of interval linear systems - overde- 
termined systems. Simply said these systems have more equations than 
C/J variables. The solution set of an interval linear system is a collection of all 

O solutions of all instances of an interval system. By the instance we mean 

a point real system that emerges when we independently choose a real 
I number from each interval coefficient of the interval system. Enclosing 

^ the solution set of these systems is in some ways more difficult than for 

QQ square systems. The main goal of this work is to present various met- 

CO hods for solving overdetermined interval linear systems. We would like to 

C~^~ present them in an understandable way even for nonspecialists in a field 

\l of linear systems. The second goal is a numerical comparison of all the 

■^h methods on random interval linear systems regarding widths of enclosures, 

C*^ computation times and other special properties of methods. 
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1 Introduction 

Real-life problems can be described by different means - by linear and nonlinear 
systems, by systems of difference and differential equations etc. Nevertheless the 
description can often be transformed to another one using only linear equalities 



(or inequalities). To take into account rounding errors or imprecisions of mea- 
surement or data we can use means of interval analysis. That is why interval 
linear systems are still in the focus of researchers. There are plenty of methods 
for enclosing the solution set of square interval linear systems - systems in the 
form Ax = b, where A is a square matrix (e.g. [3], dj, 0, [H]). That is 
because square matrices can posses some advantageous properties. They can 
be diagonally dominant, positive definite, M-matrices and many more. And we 
know that many algorithms behave well for those cases. However sometimes 
we have to work with systems called overdetermined, simply said, they consist 
of more equations than variables. That is why our favorite methods for square 
systems usually can not be applied. 

Fortunately there exist some methods for solving those systems - Gaussian 
elimination, classical iterative methods, Rohn method, least squares method or 
linear programming. The main goal of this work is to present the overview of 
existing methods for computing enclosures of solution sets of overdetermined 
interval linear systems. To the best of our knowledge, that has not been done 
for overdetermined systems yet. We would like to explain how these methods 
work in a brief but understandable way even for researchers from various fields 
not so familiar with interval linear systems. We also would like to mention some 
pitfalls and specialties connected with these methods. Some of them behave in 
an useful way - they are able to compute a very narrow enclosure or an interval 
hull or they can reveal unsolvability of an interval overdetermined system. On 
the other hand, many of them can not decide whether a solution of an interval 
system is unbounded or whether the system is unsolvable. Moreover sometimes 
if a system has certain properties, they are not able to return any meaningful 
result. It is interesting to observe how some efficient methods fail when the 
radii of intervals change and, on the other hand, to see how some simple, one 
could say "stupid" methods rule. After introducing the methods we provide 
the result of exhaustive numerical comparison of these methods concerning the 
speed, enclosure tightness and some other properties. As always it holds that 
not every method is useful for every problem case. When describing the methods 
we would like to point out the cases when they are useful. 

2 Basic notation and definitions 

It is favourable to start with the basic notation. Here we provide only a small 
part of it which we will use subsequently. For deeper introduction to interval 
arithmetics do not hesitate to see [3] . We denote interval structures in boldface 
{A,b). Point real structures are denoted in regular type {A,b). Here we work 
only with closed real intervals [c, d] , where c < d. The set M stands for the set of 
all real closed intervals. An interval x can be defined in two ways. The first one 
is by upper and lower bound x = [x, x] . The next one is by midpoint and radius 
X = {xc,x^). With intervals we can build more complex structures - interval 
matrices (of which vectors are special case). Their definition is analogous to 
the definition of a real interval {A — [A, A] or A = {Ac,Aa))- The matrix 



Ac is called the midpoint matrix and A^ is called the radius matrix. In the 
section about comparison of methods we will use the notion width of an interval 
X — [x, x] defined as w{x) = x — x. For every vector x G K" we define its sign 
vector sgn(a;) e {±1}" as 



(sgnx), = 
For a given vector x e M" we denote 



1, 
-1, 



if Xi > 0, 
if Xi < 0. 



D„ 



diag(a;i,...,a;„) 



/ xi 






X2 



^ 





V ... Xn / 

We continue with a definition of an overdetermined interval linear system. 

Definition (Overdetermined interval linear system). Let us have an in- 
terval matrix A G IM™ ", where m > n and an interval vector b G IM™. We 
call 

Ax = b 

an overdetermined interval linear system (OILS). 

Simply said an overdetermined interval linear system is just a general interval 
linear system that consist of more equations than variables. 

When we talk about the interval linear systems we have to mention what do 
we mean by the solution of an interval linear system (ILS). 

Definition (Solution set of ILS). The solution set S of an interval linear 
system Ax = b is 



T. = {x \ Ax = b fo 



A(^A,bc,b}. 



In another words it is a collection of all solutions of all instances of an interval 
linear system. By an instance we mean a point-real system that we get when we 
independently choose a real number from each interval coefficient of an interval 
system. We have to mention that this approach is different from the least 
squares approach. If no instance of the system has solution, we call the whole 
interval system unsolvable. To give a reader better conception of an interval 
solution set, on the picture IT] we show the solution set of an interval system 



-5,10]x + 10 y 
10 a; + -by 

IQx + [10,25]j/ 



[15,20]z = [50,100], 

[5,15]z = [-50,50], 

[-10, -5] z = [50,100]. 



The solution set is generally a polyhedral set, not necessarily convex. Never- 
theless it is convex in each orthant. As we can see, this set is quite difficult to 




Figure 1: Graphical view of a solution of an interval linear system with 3 vari- 
ables 



be described. There exist a lot of possibilities of description. One of them is to 
find an n-dimensional box (aligned with axes) as tight as possible that contains 
the solution set. We call this box the interval hull. Unfortunately it can be 
proved that computing the interval hull is NP-hard. That is why we look for 
a box a little bit wider, still as narrow as possible, that contains the interval 
hull. We call this box the interval enclosure. There is a large variety of methods 
we can use to compute an interval enclosure. We have already built sufficient 
terminology, therefore we are able to move on to description of methods. 

3 Gaussian elimination (GE) 

The first thing that can come to our mind is the Gaussian elimination. The 
interval GE was proposed by Hansen in [1]. The idea is pretty the same as for 
point real systems. We present a slightly modified version. We eliminate rows 
using interval operations with the only difference that we eliminate the system 
{A\b) to the shape 

where C is an (n — 1) x (n — 1) interval matrix in the Row echelon form (REF) 
with [1, 1] intervals on pivot positions, d, e are (n — 1) x I interval vectors, is 
an {m — n + 1) x [m — n + 1) matrix composed of the intervals [0, 0] and /,g 
are {m — n + 1) x 1 interval vectors. 

One thing is often not clear to people familiar with interval arithmetics when 
Gaussian elimination is explained to them. And that is why on pivot positions 
there are [1,1] intervals and why there are [0,0] intervals in matrix. This 
might not be completely obvious. However the reason is quite simple - we can 
think of working with pivots and elements beneath as performing the elimina- 
ting operations on them in all instances of an interval system separately (after 
elimination pivots in instances are equal to 1, elements under the pivots are 



equal to 0). When we assemble the elimination result, we know that we get 
[1, 1], [0,0] intervals on these positions in all these instances, therefore there is 
no need to overestimate them using interval arithmetics. For the other coeffi- 
cients we have to use interval arithmetics to gain their interval enclosures. 

Now we can realize that vectors /,g form m — n + I interval equations in the 
shape 

f,x„ = g^ for i = 1, . . . , (m - n + 1). 

The solution of these equations is intersection of all intervals Qi/ fi - we get an 
enclosure of the variable a;„. If the intersection is empty, then the system has no 
solution. Nonetheless, if the intersection is unbounded, it can either mean that 
the solution set of the system is unbounded or large overestimation due to large 
number of interval operations occurred. The enclosures for the other variables 
can be obtained using backward substitution as we know it (only computing in 
interval arithmetics). 

This algorithm works only for very small m x n systems (n ~ 4). For larger 
systems we get a huge overestimation as we can see in table 111 



variable 


midpoint 


radius 


Xi 


0.1479 


227.6698 


X2 


15.2091 


172.5929 


X3 


11.1031 


68.4653 


X4 


9.7809 


64.1056 


Xr, 


-8.8168 


27.2234 


X6 


25.8164 


25.8398 


X-j 


-19.0444 


30.4596 


Xs 


-22.0799 


11.0313 


Xg 


1.9649 


12.1172 


a^io 


-19.1817 


11.6841 


Xu 


-20.9670 


1.9153 


X\2 


-4.6988 


3.5407 


X\z 


3.1223 


4.3894 



Table 1: The overestimation of the enclosure (random system 15 x 13) with 
random radii < 10"'^ caused by GE and backward substitution 

For larger systems it is necessary to use preconditioning. We use the pre- 
conditioner mentioned by Hansen in [T]. 



C 



A\ 
A% I 



where A^ consists of first n rows of Ac (approximation of a midpoint matrix 
of A) and A2 consists of the remaining m — n rows of Ac. Operation ~^ is an 
approximate inverse operation. Of course a preconditioning usually causes an 
overestimation of resulting enclosure, but it often causes smaller enlargement 



compared to performing GE without preconditioning. We must also mention 
that preconditioning enlarges the solution set of the system, that is why un- 
solvable systems become solvable after preconditioning. Therefore if we want to 
check unsolvability we can not use preconditioning. If we use the elimination to 
Reduced row echelon form (RREF) and allow multiple right-hand side, we are 
able to compute an interval matrix containing A~^. For a proper definition of 
the interval inverse matrix see [7]. 

4 Iterative methods 

There are many interval modifications of methods for square non-interval sys- 
tems - Jacobi, Gauss-Seidel, Krawczyk, Hull method. But we usually can not 
use them when dealing with overdetermined systems. Moreover without precon- 
ditioning these methods often work really badly. However we have one possibil- 
ity, if we want to use our known iterative methods we can use the preconditioning 
mentioned in the section about Gaussian elimination. We can see that after this 
preconditioning the new system is usually of the shape 

0" 

where J~ is an n x n interval matrix "wrapping" the point eye matrix and 0~ 
is an {m — n) x n interval matrix with small intervals containing zero. Then we 
can use our favorite iterative method for the upper square nx n subsystem. Let 
us choose Jacobi method for example. We start with some initial enclosure a;(°) 
and iteratively "sharpen" this enclosure using the following formula (in each 
step we simply express the variable Xi from i-th equation, for more information 
see [3]) 

X* = -^{b^ - Y^ A^jxf^) for i = (1, . . . ,n), in (fc + l)-th step. (1) 
After every iteration we provide intersection with the old enclosure to obtain 

(fe+l) (fc) _ * r . /, X 

x] ' — XI ' r\x^ for I = (1, ... , n). 

Gauss-Seidel method differs in only one thing - we do not wait and compute 
intersection after every computation of x* and immediately use the new value 
for further computing. This leads to less iteration steps. 

The problem is that when using these methods, intervals containing zero are 
are not allowed on a diagonal (since there is the division with An). That 
case happens if the radii of the original matrix are "large" (even 10~^). This 
method looks rather simple. We loose some information due to precoditioning 
and chopping the last m — n rows of the original system. But we can realize 



that e.g. Jacobi method can be effectively parallehzed. And if there is effective 
way to determine x^^\ then this method can be used for very fast sharpening of 
x^^> , for example in Constraint Programming as a sharpening step in between 
the run of another sharpening methods. 

5 Rohn method 

We would like to mention the method introduced by Rohn. Due to lack of space 
we will introduce it only briefly. For more information and theoretical insight 
one can take a look in [S] . The base of the method lies on the following theorem. 

Theorem (Rohn). Let Ax ^ b be an lOLS with a solution set E (A is an 
m X n matrix). Let R be an arbitrary real n x m matrix, let xq and d > be 
arbitrary n-dimensional real vectors such that 

Gd + g < d, 

where 

G = \L - RAc\ + \R\Aa 

and 

g = \R{AcXo - 6c)| + \R\{Aa\xo\ + 6a). 

Then 

S C [a:Q — d, Xq + d]. 

The question is how to find the vector d, the matrix R and the vector Xo- To 
compute d, we can, for example, rewrite the inequality as 

d = Gd + g + f, 

for some small vector / > 0. Then start with d ~ and iteratively refine d. 
This algorithm will stop after a finite number of steps if the spectral radius of 
G is less than 1. Otherwise we do not know what happens. During practical 
testing with random systems with radii of intervals close to 0.1 the vector d > 
was rarely found. 

We still have to determine xq and R. For the start we can take 

Xq « Rbc, 

R^{AlA^)-'Al 

but not necessarily. The theorem provides a clever instrument for iterative 
improvement of solution. We do not have to use only Ac to compute R, we can 
take any (e.g. random) A € A, compute an enclosure and then intersect it with 
the old one. We can repeat this process as many times we want and provide 
an iterative improvement of the enclosure. For smaller systems the iterative 



improvement works pretty well. The narrowing of initial enclosures is showed 
in the table [2J It shows the ratios of enclosure widths returned by non-iterative 
Rohn method and iterative versions for 10, 100, 1000 iterations. To see how 
the ratios are computed, it is possible to take a look in subsection 8.1. In the 
further text we will call the first method basic Rohn method and use the notion 
basic enclosure. When we use the iterative improvement we will talk about the 
iterative Rohn method and iterative enclosure. 



system 


10 it. 


100 it. 


1000 it. 


5x3 


0.73 


0.57 


0.50 


15 X 10 


0.89 


0.82 


0.76 


25 X 21 


0.94 


0.90 


0.87 


35 X 23 


0.95 


0.92 


0.90 


50 X 35 


0.97 


0.94 


0.93 


70 X 55 


0.98 


0.96 


0.95 


100 X 87 


0.98 


0.97 


0.97 



Table 2: Rohn method - ratios of widths of basic enclosures and iterative en- 
closures (10, 100 a 1000 iterations) 



6 Least squares method 

We can also make use of the least squares formula A^ Ax = A^b. But when we 
write it in an interval way 

A^Ax = A^b 

that indeed does not work because of two interval matrix multiplication. Even 
if we use some preconditioner C and write 

{CAf{CA)x = {CAfb 

that does not work either. Anyway we can use an equivalent expression of the 
least squares formula for point real systems. If we use it for interval systems we 
get 

/ A\ f y \ f b 











Now we have a square system and we can apply some suitable method for square 



systems. We get the vector {y,x)^ G 



j)m+n 



as the result and take x as the 



solution enclosure. It can be seen that the returned interval vector contains the 
solution of the interval least squares, therefore this method returns a solution 
even if the system is unsolvable. Another drawback is that if the original system 
is of size to x n we have to solve a new one of size (m + n) x (m + n) On the 
other hand this method computes very sharp interval enclosures. It is used in 
the method verifylssO in INTLAB 6 [TT] in combination with the interval 



Krawczyk method (described e.g. in [3]) for square systems. When solving a 
systems in the shown way, we can notice the dependencies in the new system. We 
can possibly use methods dealing with dependencies in interval linear systems 

(e.g. m, El, m)- 

7 Linear programming (LP) 

It is possible to compute the interval hull by using linear programming. We 
can use the famous theorem by Oettli-Prager. It can be found with proof for 
example in [8]. Due to lack of space, here we omit the proof. 

Theorem (Oettli-Prager). Let us have an interval linear system Ax = b. 
Vector X £ M" is a solution of this system (x € T,) if and only if 

\AcX-bc\ < AaI^I +6a- 

Unfortunately we are still not able to use LP because of the absolute values. 
Now we demonstrate how to rewrite this problem using linear inequalities only. 
We can get rid of the first one by decomposing it in two cases 

AcX-bc <Aa\x\ +6a, (2) 

-(AcX-bc) <AA\x\+bA. (3) 

The second absolute value can be rewritten with the use of knowledge of the 
orthant we currently "are" inside. The following holds 

\x\ — DzX, where z — sgn(a;). 

That gives a rise to the condition 

< D,x. (4) 

Now we can see, that for every orthant conditions (pi), (Is]) and Q form a system 
of linear inequalities. Therefore we can use linear programming. Unfortunately 
we have to use (2" x 2n) linear programming problems (in each coordinate we 
compute the upper and lower bound). That is obviously too much computing. 
However we can compute the enclosure of the solution of the system with some 
suitable method (least squares, Rohn) and than apply linear programming to 
only those orthants where this enclosure lies. This approach is often much 
more faster. The table [3] illustrates a large speedup when solving the system 
with different method before applying LP. The sign '-' means that we omitted 
the testing because of enormous computing time when compared to LP with 
presolving. 

8 Comparison of methods 

The another goal of this paper is a numerical comparison of methods for enclos- 
ing solutions of overdetermined systems. As far as we are concerned for these 



system 


time LP 


time LP presolved 


5x3 


6 sec 


1 sec 


9x5 


43 sec 


1.68 sec 


13 X 7 


5 min 


3.59 sec 


15 X 9 


28 min 


4.1 sec 


25 X 21 


- 


13 sec 


35 X 23 


- 


19 sec 


45 X 31 


- 


43 sec 


55 X 35 


- 


1 min 


73x55 


- 


9 min 



Table 3: Comparison of times of LP in all orthants and LP with presolving with 
a different method (here verifylss) 

methods it has not been done yet. All subsequent tests were computed using 
the following hardware and software: 

• Processor - AMD Phcnom(tni) II X6 1090T 

• Memory - 15579 MB 

• Matlab R2010b 

• Intlab 6 (see [llj ) 

• Versoft 10 (see [TD]) for linear programming 

It is difficult to imagine rectangular matrices in some special shape, therefore 
we test the methods on random matrices. In this work we use four parameters: 

• maximum radius 

• m.idpoint range 

• stopping parameter 

• maximum number of iterations 

Most of these parameters are used for generating random interval systems 
with certain properties. The key parameter is maximum radius of interval co- 
efficients of an interval linear system. When this parameter is for example 
—4, that means that for every interval coefficient y = (j/c,2/a) in the system 
j/A £ 10"'' holds. It is usually negative, because radii greater than 1 will lead to 
singularities. The radii in a system are not the same, they are chosen randomly 
with respect to the maximum radius parameter. Another important parameter 
is the midpoint range, which specifies the range of midpoints of intervals. The 
methods may have different properties when working with all intervals that have 
relatively close midpoints and between the case when the midpoints are quite 
far like in vector (10, 1, 11234, 0.1)"^. We will test on two cases - the first one is 



10 



with midpoints uniformly chosen from [—25, 25] and the second one uniformly 
chosen from [—1000, 1000]. Some methods need small real positive number e as 
a stopping parameter (some methods need a small positive vector, but we can 
use / = (e, e, . . . , e)^). We will choose e = lo(™=^=^™"'" radius-2)_ ^f^g^. perform- 
ing some practical tests this seems to be a reasonable choice. Some iterative 
methods use parameter maximum of iterations. It is a safety precaution for 
the case the stopping parameter did not work. We use the value 20 for all the 
tested cases. We will test systems up to size about m = 1000 which equals to 
a thousand of linear equations. We chose this upper bound on size of tested 
systems after consulting with a colleague from the Faculty of Civil Engineering. 
These are maybe the largest linear systems they need to solve for many technical 
purposes. If a solution of larger system is needed, then it is usually better to 
split the problem into parts, solve the parts separately and then assemble the 
solution. For the sake of clarity we establish the following identifiers of methods: 

• GSpre - Gauss-Seidel iterative method with preconditioning 
GEpre - Gaussian elimination with preconditioning 

• Rohn - basic Rohn method 
Verifylss - method vcrifylss from toolbox Intlab by Rump 

• LPver - linear programming from toolbox Versoft by Rohn 

8.1 Enclosure width comparison 

We are going to compare the widths of enclosures of all the mentioned meth- 
ods according to test enclosures returned by a different method. For an n- 
dimensional enclosure vector x = (a;i, 0:2, . . . , Xn)^ and a;**^"* as a test enclosure 
we compute the ratio 



• 



• 



ratio(a;,a;*''^*) 



1 .^ w(a;^ 



, w(a;* 



Clearly, the lower the ratio, the better the enclosure. 

We divide the testing in two phases. For larger systems it is very time 
consuming to compute many enclosures using LPver, moreover for small radii 
of intervals the method returns NaN solution. That is why we tested the 
methods on small systems first and the enclosures were compared with LPver 
results and in second phase they were compared to Verifylss enclosures. First 
let us show the table |4] where the ratios of small systems are displayed. 

We can see that GEpre and GSpre are doing really badly. And that their 
enclosure ratio is of about the same value. We attribute that to the use of the 
same preconditioner. Rohn and Verifyllss provide good results. Verifylss is 
a little bit better, but the difference is slowly disappearing with the grow of a 
system size. That is confirmed by the second table[5] We made the radii smaller 
(< 10^^) since GEpre often returned infinite solution due to a large number of 
interval operations. 
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system 


GSpre 


GEpre 


Rohn 


Verifylss 


5x3 


2.9968 


2.9969 


1.2347 


1.1893 


15x9 


7.7057 


7.7069 


1.1601 


1.1500 


35 X 23 


8.3591 


8.3602 


1.1276 


1.1249 


55 X 35 


11.5026 


11.5064 


1.1336 


1.1331 


73 X 55 


17.6528 


17.6624 


1.0828 


1.0848 



Table 4: Ratios of widths of enclosures, LPver is 1 (maximum radii < 10 ) 



system 


GSpre 


GEpre 


Rohn 


100 X 45 


21.1668 


21.1668 


1.0237 


100 X 87 


9.2420 


9.2420 


1.0062 


180 X 125 


26.6744 


26.6744 


1.0064 


180 X 170 


15.0262 


15.0262 


1.0020 


290 X 190 


37.6976 


37.6977 


1.0044 


290 X 260 


26.9361 


26.9361 


1.0018 


380 X 275 


40.5159 


40.5160 


1.0028 


380 X 360 


32.9500 


32.9504 


1.0009 


500 X 350 


66.9764 


66.9771 


1.0022 


500 X 470 


34.2512 


34.2516 


1.0007 



Table 5: Ratios of widths of enclosures, Verifylss is 1 (maximum radii < 10 



8.2 Enclosure computation time comparison 

We now compare the computation times of methods - GSpre, GEpre, Rohn, 
Verifylss. Midpoint matrices and vectors are chosen uniformly from the interval 
[-1000, 1000]. We provide two tables p^ and ^ for maximum radius < 10~^ and 
< 10~^ respectively. They display the average computation times of all methods 
for systems of different sizes. Time is measured with Matlab functions tic and 
toe. The results are displayed only to the fourth decimal number. More digits 
are beyond the precision of tic and toe. 



system 


GSpre 


GEpre 


Rohn 


Verifylss 


5x3 


0.0113 


0.0117 


0.00065 


0.0021 


15 X 13 


0.0367 


0.0685 


0.00086 


0.0023 


35 X 23 


0.0636 


0.2840 


0.0012 


0.0036 


50 X 35 


0.0983 


0.5806 


0.0019 


0.0052 


100 X 87 


0.3063 


2.4631 


0.0115 


0.0169 


200 X 170 


0.6632 


9.8570 


0.0523 


0.0761 


380 X 275 


1.3253 


34.5861 


0.1765 


0.3257 


500 X 470 


2.9486 


69.6341 


0.6235 


0.9154 



Table 6: Average times (in seconds) of all methods for maximum radius < 10 '^ 
We can see that the methods do not show any remarkable sensibility to the 
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system 


GSpre 


GEpre 


Rohn 


Verifylss 


5x3 


0.0113 


0.0118 


0.00066 


0.0021 


15 X 13 


0.0356 


0.0687 


0.00086 


0.0023 


35 X 23 


0.0608 


0.2836 


0.0012 


0.0036 


50 X 35 


0.0911 


0.5809 


0.0019 


0.0053 


100 X 87 


0.2416 


2.5224 


0.0117 


0.0171 


200 X 170 


0.5718 


10.6842 


0.0561 


0.0760 


380 X 275 


1.2515 


37.0210 


0.1924 


0.3427 


500 X 470 


2.5774 


75.5171 


0.6616 


0.9761 



Table 7: Average times in seconds of all methods for maximum radius < 10 



-5 



changes of radii of intervals if we keep the stopping criterion equal to (maximum 
radius —2). There was also no extraordinary change when lowering the midpoint 
range to something much smaller (e.g. [—25,25]). The complete time winner 
is Rohn, the second is Verifylss. The method GEpre is doing really badly. 
We shall look more carefully on the comparison between Rohn and Verifylss. 
The method Rohn is almost always faster. According to the table [5] we can 
object that it pays off to take a little bit more time and compute more exact 
enclosure using Verifylss. But let us look closer on the following table |8] It 
displays the ratios of computation times of Rohn and Verifylss for systems of 
broader spectrum of sizes. We did the tests for many random systems, here we 
selected some interesting ones. It is possible to see, that for m x n systems, 
where m is much bigger than n, Rohn is much more faster. For example on 
system 278 x 35 we can call the iterative Rohn with 13 iterations within the 
same time Verifylss needs to do the computation. And we can possibly sharpen 
the enclosure. The question whether this procedure pays off is answered by the 
table [2J For relatively small systems it definitely pays off. For larger systems 
we might get not so significant improvements and it is safer to use Verifylss. 



system 


t{Verifylss) 
t(Rohn) 


80 X 41 


2.8 


189 X 166 


1.5 


278 X 35 


13.8 


377 X 319 


1.6 


525 X 285 


2.7 


712 X 271 


4.5 


807 X 68 


46.5 


894 X 8 


423.7 


894 X 797 


1.5 


906 X 128 


19.1 


978 X 235 


9.4 


1000 X 663 


2.0 



Table 8: Ratios of Verifylss and Rohn times 
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9 Conclusion 

We introduced many methods for solving overdetermined interval linear systems 
- Gaussian elimination, Iterative methods, Rohn method. Least square method 
and Linear programming. We mentioned some advantageous and unfavourable 
properties of these methods. At the end of this work we showed the comparison 
of all mentioned methods. As usually it is hard to tell which method is the 
best. Linear programming is suitable when we want to compute the interval 
hull of an interval system and have enough time. When our system has really 
small coefficients or our system is in the shape of " noodle" , it is advantageous 
to use Rohn method otherwise it is favourable to use the least square method 
for fast computing of a sharp interval enclosure. Iterative methods can be used 
when we already know some enclosure and want to make try to sharpen it (e.g. 
as a part of some constraint solver) and then passing it as an input to another 
method c.g LP. Gaussian elimination seems to work really bad, but when we 
consider the solvability of the system, it might give us a valuable information 
about the nonexistence of solution. 

Computing enclosures of solutions of overdetermined interval linear systems 
can occur as a subproblem of many computational tasks (constraint satisfaction 
problem). And also the information about solvability of a system is often im- 
portant (system validation, technical computing). There is still much work to 
do in this area. We believe that even faster and sharper algorithms can be de- 
veloped. Simultaneously we are working on methods checking the solvability of 
overdetermined interval linear systems. Hopefully sufficient or necessary condi- 
tions for solvability or unsolvability of OILS can be found. The preconditioning 
is used for many methods, it would be useful to test different types of precon- 
ditioners. Some methods use probabilities, e.g. iterative Rohn method. There 
is a question whether we can effectively derandomize these methods. 
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